Stretching dependence of the vibration modes of a single-molecule Pt-H 2 -Pt bridge 
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A conducting bridge of a single hydrogen molecule between Pt electrodes is formed in a break 
junction experiment. It has a conductance near the quantum unit, Go = 2e 2 /h, carried by a 
single channel. Using point contact spectroscopy three vibration modes are observed and their 
variation upon isotope substitution is obtained. The stretching dependence for each of the modes 
allows uniquely classifying them as longitudinal or transversal modes. The interpretation of the 
experiment in terms of a Pt-Eh-Pt bridge is verified by Density Functional Theory calculations for 
the stability, vibrational modes, and conductance of the structure. 

PACS numbers: 73.63.Rt, 63.22. +m, 73.23.-b, 85.65.+h 



There is beauty and power in the idea of constructing 
electronic devices using individual organic molecules as 
active elements. Although the concept was proposed as 
early as 1974 Q only recently experiments aimed at con- 
tacting individual organic molecules are being reported 
|J|JlllllIllEllIl3,II3 and devices are being tested 
12l|lj|. The first results raised high expectations, but 
quickly problems showed up such as large discrepancies 
between the current-voltage characteristics obtained by 
different experimental groups, and large discrepancies be- 
tween experiments and theory. The main tools that have 
been applied in contacting single molecules are STM (or 
conducting tip AFM) and break junction devices. Of- 
ten it is difficult to show that the characteristics are 
due to the presence of a molecule, or that only a sin- 
gle molecule has been contacted. There has been im- 
portant progress in analysis and reproducibility of some 
experiments 0, IE IS El> Ell 0- but in comparing the 
data with theory many uncertainties remain regarding 
the configuration of the organic molecule and the nature 
of the molecule-metal interface. The organic molecules 
selected for these studies are usually composed of several 
carbo-hydride rings and are anchored to gold metal leads 
by sulphur groups. In view of the difficulties connected 
with these larger molecules it seems natural to step back 
and focus on even simpler systems. 

Here we concentrate on the simplest molecule, H2, an- 
chored between platinum metal leads using mechanically 
controllable break junctions. The first experiments on 
this system [l4| showed that the conductance of a sin- 
gle hydrogen molecule between Pt leads is slightly below 
IGo, where Go = 2e 2 /h is the conductance quantum. 
A vibration mode near 65 meV was observed and in- 
terpreted as the longitudinal center-of-mass (CM) mode 
of the molecule. These results have inspired new calcu- 
lations on this problem using Density Functional The- 
ory (DFT) methods [3 El. Cuevas et al. find a 
conductance around 0.9Go, in agreement with the first 



DFT calculations presented in In contrast, Garcia 
et al. 16] obtain a conductance of only (0.2 — 0.5) Go 
for the in-line configuration of the hydrogen molecule. 
Instead they propose an alternative configuration with 
hydrogen atoms sitting above and below a Pt-Pt atomic 
contact. 

In this letter we combine new experimental results with 
DFT calculations to show that the configuration pro- 
posed in E3 is correct, yet the observed vibration mode 
was incorrectly attributed. In contrast, the present ex- 
periment resolves three vibration modes that can be clas- 
sified as longitudinal or transverse modes based on the 
observed shifts with stretching of the contacts. The com- 
parison with the calculations is nearly quantitative and 
the large number of experimentally observed parameters 
(the number of vibration modes, their stretching depen- 
dence and isotope shifts, the conductance and the num- 
ber of conductance channels) puts stringent constraints 
on any possible interpretation. 

The measurements have been performed using the me- 
chanically controllable break junction technique [l7lll8|. 
A small notch is cut at the middle of a Pt wire to fix the 
breaking point. The wires used are 100 /^m in diameter, 
about 1 cm long, and have a purity of 99,9999%. The 
wire is glued on top of a bending beam and mounted in a 
three-point bending configuration inside a vacuum cham- 
ber. Once under vacuum and cooled to 4.2 K the wire 
is broken by mechanical bending of the substrate. Clean 
fracture surfaces are exposed and remain clean for days 
in the cryogenic vacuum. The bending can be relaxed to 
form atomic-sized contacts between the wire ends using 
a piezo element for fine adjustment. 

After admitting a small amount (~ 3 /xmole) of molec- 
ular H2 (99.999%) in the sample chamber and waiting 
some time for the gas to diffuse to the cold end of the 
insert, a sudden change is observed in the conductance 
of the last contact before breaking. The typical value of 
(1.6 ± 0.4)Gq for a single-atom Pt contact is replaced by 
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FIG. 1: Differential conductance curve for D2 contacted by 
Pt leads. The dl/dV curve (top) was recorded over 1 minute, 
using a standard lock-in technique with a voltage bias mod- 
ulation of 1 meV at a frequency of 700Hz. The lower curve 
shows the numerically obtained derivative. The spectrum for 
H2 in the inset shows two phonon energies, at 48 and 62 meV. 
All spectra show some, usually weak, anomalies near zero bias 
that can be partly due to excitation of modes in the Pt leads, 
partly due to two- level systems near the contact JT3| . 



tance, in which case we observe an increase of the vibra- 
tion mode energy with stretching. This observation suf- 
fices to invalidate the original interpretation of this 
mode as the longitudinal CM mode. Indeed, our DFT 
calculations show that the stretching mainly affects the 
Pt-H bond which is elongated and weakened resulting 
in a drop in the frequency of the H2 longitudinal CM 
mode. An increase can be obtained only for a transverse 
mode which, like a guitar string, obtains a higher pitch at 
higher string tension due to the increased restoring force. 

On many occasions we observe two modes in the dl / dV 
spectra, see the inset of Fig. ^ The relative amplitude 
of the two modes varies; some spectra show only the 
lower mode, some only the higher one. All frequencies 
observed in a large number of experiments are collected 
in the histograms shown in Fig. [21 With a much larger 
data set compared to Q we are now able to resolve 
two peaks in the distribution for H2 corresponding to 
the two modes seen in the inset of Fig. ^ The peaks 
are expected to shift with the mass m of the isotopes as 
uj oc y/l/m. This agrees with the observations, as shown 
by the scaled position of the hydrogen peaks marked by 
arrows above the distributions for D2 and HD. Note that 
the distribution for HD proves that the vibration modes 
belong to a molecule and not an atom, since the latter 



a frequently observed plateau near IGo that has been at- 
tributed in [l4| to the formation of a Pt-H2-Pt bridge. By 
increasing the bias voltage above 300 mV we recover the 
pure Pt conductance. But as soon as the bias voltage 
is decreased the H2 induced plateaus at IGo reappear. 
We interpret this as desorption of hydrogen due to Joule 
heating of the contacts. For biases below 100 mV, the 
Pt-H2-Pt bridge can be stable for hours. 

At the IGo-conductance plateaus we take differential 
conductance (dl/dV) spectra in order to determine the 
inelastic scattering energies. By repeatedly breaking the 
contacts, joining them again to a large contact, and 
pulling until arriving at a plateau near IGo, we obtain 
a large data set for many independent contacts. The 
experiments were repeated for more than 5 independent 
experimental runs, and for the isotopes HD (96%) and 
D2 (99.7%). Fig.f^shows a spectrum taken for D2 show- 
ing a sharp drop in the differential conductance by 1-2% 
symmetrically at ±50 meV. Such signals are character- 
istic for point contact spectroscopy I20J , which was first 
applied to single-atom contacts by |2lJ . The principle of 
this spectroscopy is simple: when the difference in chem- 
ical potential between left- and right-moving electrons, 
eV, exceeds the energy of a vibration mode, tkJ, back- 
scattering associated with the emission of a vibration be- 
comes possible, giving rise to a drop in the conductance. 
This can be seen as a dip (peak) in the second derivative 
d 2 I/dV 2 at positive (negative) voltages, as in Fig.fl] 

Some contacts can be stretched over a considerable dis- 
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FIG. 2: Distribution of vibration mode energies observed for 
H2, HD, and D2 between Pt electrodes, with a bin size of 
2 meV. The peaks in the distribution for H2 are marked by 
arrows and their widths by error margins. These positions 
and widths were scaled by the expected isotope shifts, \/2/3 
for HD and \/T/2 for D2, from which the arrows and margins 
in the upper two panels have been obtained. 
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would have produced a mixture of the H2 and D2 distri- 
butions. In the case of D2 we observe a third peak in the 
distribution at 86-92 meV. For the other isotopes this 
mode falls outside our experimentally accessible window 
of about ±100 meV, above which the contacts are desta- 
bilized by the large current. Fig. shows the dependence 
of this mode upon stretching of the junction. In contrast 
to the two low-frequency modes this mode shifts down 
with stretching, suggesting that this could be the longi- 
tudinal CM mode that was previously attributed to the 
low- frequency modes |l4|. 

In order to test the interpretation of the experiment in 
terms of a Pt-tb-Pt bridge we have performed extensive 
DFT calculations using t he p lane wave based pseudopo- 
tential code Dacapo |22j, |2j| . The molecular contact is 
described in a supercell containing the hydrogen atoms 
and two 4-atom Pt pyramids attached to a Pt(lll) slab 
containing four atomic layers, see inset of Fig. 0] [2J] . In 
the total energy calculations both the hydrogen atoms 
and the Pt pyramids were relaxed while the remaining 
Pt atoms were held fixed in the bulk structure. The vi- 
bration frequencies were obtained by diagonalizing the 
Hessian matrix for the two hydrogen atoms. The Hes- 
sian matrix is defined by d 2 Eq / {du na du rn p) , where Eq is 
the ground state potential energy surface and u na is the 
displacement of atom n in direction a multiplied by the 
mass factor \JM n . In calculating the vibration modes all 
Pt atoms were kept fixed which is justified by the large 
difference in mass between H and Pt. The conductance 
is calculated from the Meir-Wingreen formul a 1271 using 
a basis of partly occupied Wannier functions |28j , repre- 
senting the leads as bulk Pt crystals. 
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FIG. 3: dl /dV Spectrum (top) for Pt-D2 junctions, and their 
numerical derivatives d 2 I/dV 2 (bottom). The spectrum (b) 
was obtained on the same contact as (a) after stretching the 
junction by 0.04 nm. They are two spectra out of a series of 
six and the complete evolution of the two modes is shown in 
the inset. 




FIG. 4: Calculated vibrational frequencies for the hydrogen 
atoms in the contact as a function of elongation of the super- 
cell. The inset shows the atomic arrangement in the supercell. 
The lower panel shows the calculated conductance. 

In order to simulate the stretching process of the ex- 
periment we have calculated contacts for various lengths 
of the supercell. The bridge configuration is stable over 
a large distance range with the binding energy of the H2 
molecule varying from —0.92 eV to —0.47 eV, relative to 
gas phase H 2 and a broken Pt contact, over the range of 
stretching considered here. The H-H bond length stays 
close to 1.0 A during the first stages of the stretching 
upon which it retracts and approaches the value of the 
free molecule. The hydrogen thus retains its molecular 
form and the elongation mainly affects the weaker Pt- 
H bond. For smaller electrode separations a structure 
with two hydrogen atoms adsorbed on the side of a Pt- 
Pt atomic contact becomes the preferred geometry, as 
also found by Garcia et al. 16]. However, we find that 
the latter structure has a conductance of 1.5 Go, well 
above 1 Gq. Moreover, this structure has at least three 
conduction channels with significant transmission, which 
excludes it as a candidate structure based on the analysis 
of conductance fluctuations in an d HU , which find a 
single channel only. In view of the activity of the Pt sur- 
face towards catalyzing hydrogen dissociation one would 
have anticipated a preference for junctions of hydrogen in 
atomic form. However, we find that the bonding energy 
of H compared to that of H2 strongly depends on the 
metal coordination number of the Pt atom. For metal 
coordination numbers smaller than 7 bonding to molec- 
ular hydrogen is favored, the bond being strongest for 
fivefold coordinated Pt. 

The calculations identify the six vibrational modes of 
the hydrogen molecule. For moderate stretching two 
pairwise degenerate modes are lowest in frequency The 
lowest one corresponds to translation of the molecule 
transverse to the transport direction while the other one 
is a hindered rotation mode. The two modes are charac- 
terized by increasing frequencies as a function of stretch 
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of the contact. At higher energies we find the two lon- 
gitudinal modes: first the CM mode and then the inter- 
nal vibration of the molecule. These two modes become 
softer during stretching up to Pt-H bond lengths of about 
1.9 A. Beyond this point the Pt-H bond begins to break 
and the internal vibration mode approaches the one of 
the free molecule. 

The variation of the frequencies of the lowest lying hy- 
drogen modes with stretching is thus in qualitative agree- 
ment with the experiments, a strong indication that the 
suggested structure is indeed correct. The agreement is 
even semiquantitative: If we focus on displacements in 
the range 1.7 - 2.0 A (see Fig. 4) the calculated conduc- 
tance does not deviate significantly from the experimen- 
tally determined value close to 1G - In this regime the 
three lowest calculated frequencies are in the range 30- 
42 meV, 64-92 meV, and 123-169 meV. The two lowest 
modes can be directly compared with the experimentally 
determined peaks at 54 and 72 meV observed for H 2 , 
while a mass re-scaling of the D 2 result for the highest 
mode gives approximately 126 meV. 

The second peak in the HD distribution in Fig. [21 is 
somewhat above the position obtained by scaling the H2 
peak by \/3/2. The transverse translation mode and 
the hindered rotation mode are decoupled when the two 
atoms of the molecule have the same mass. In the case 
of HD they couple and the simple factor does not hold. 
Having identified the character of these modes a proper 
re-scaling of the experimentally determined H 2 frequen- 
cies (54 and 72 meV) to the case of HD leads to the 
frequencies 42 and 66 meV, in very good agreement with 
the peaks observed for HD. 

Even though there is good agreement with the cal- 
culated signs of the frequency shifts with stretching for 
the various modes, there is a clear discrepancy in magni- 
tudes. Considering, e.g., the high- lying mode for D 2 , the 
measured shift of the mode is of the order 15 me v/A, 
which is almost an order of magnitude smaller than the 
calculated variation of around 130 meV/A. However, ex- 
perimentally the distances are controlled quite far from 
the molecular junction and the elastic response of the 
electrode regions have to be taken into account. Simula- 
tions of atomic chain formation in gold |30| during con- 
tact breaking show that most of the deformation happens 
not in the atomic chains but in the nearby electrodes. A 
similar effect for the Pt-H 2 -Pt system will significantly 
reduce the stretch of the molecular bridge compared to 
the displacement of the macroscopic electrodes. 

The observation of the three vibration modes and their 
stretching dependence provides a solid basis for the in- 
terpretation. The fourth mode, the internal vibration, 
could possibly be observed using the isotope tritium. The 
hydrogen molecule junction can serve as a benchmark 
system for molecular electronics calculations. The ex- 
periments should be gradually expanded towards more 
complicated systems and we have already obtained pre- 



liminary results for CO and C 2 H 2 between Pt leads. 
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